Universal Distribution of Kondo Temperatures in Dirty Metals 
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Kondo screening of diluted magnetic impurities in a disordered host is studied analytically and 
numerically in one, two and three dimensions. It is shown that in the T K — > limit the distribution 
of Kondo temperatures has a universal form, P(T K ) ~ T^ a that holds in the insulating phase and 
persists in the metallic phase close to the metal insulator transition. Moreover, the exponent a 
depends only on the dimensionality. The most important consequence of this result is that the T 
dependence of thermodynamic properties is smooth across the metal-insulator transition in three 
dimensional systems. 

PACS numbers: 71.27. +a, 75.20.Hr, 71.10.Hf 



In a class of compounds with coexisting magnetic mo- 
ments and conduction electrons the magnetic susceptibility 
and the specific heat coefficient exhibit anomalous power 
law or logarithmic low temperature divergences. However, 
magnetic order is not observed down to the lowest attain- 
able temperatures . This behavior is observed in a grow- 
ing family of heavy-fermion compounds Q an also in the 
doped semiconductor Si:P 0. 

It was early recognized that this type of non-Fermi liq- 
uid (NFL) anomaly is associated with the interplay of the 
Kondo effect and strong disorder. In heavy fermions dis- 
order is present mostly in the /-electron sub-system. In a 
clean system screening of /-moments starts at the single- 
impurity Kondo scale T K ~ exp[— 1/( Jpo)] where J is the 
Kondo coupling and po is conduction electron density of 
states (DOS) at the Fermi level. In the dirty case J is a 
random variable with a distribution P(J). Therefore, the 
local values of T K are also random and have a distribu- 
tion P(T K ) solely determined by P{J). The simple Kondo- 
disorder model 0] and its subsequent improvements []|, 
based on the above physical picture, lead quite naturally 
to NFL phenomenology. 

One might think that when the disorder is mostly present 
in the conduction band P(T K ) is similarly determined by 
p(R,e F ), the distribution of the local DOS (LDOS) of the 
host at the Fermi level at the position of the Kondo impu- 
rity 0. 

In this Letter we show that this is not the case except 
at weak disorder. In general, P(T K ) depends on subtle fea- 
tures of the LDOS and its statistical properties, not sim- 
ply on the distribution of its values at the Fermi energy 
13, la, |9j . The most important consequence of the analysis 
presented below is that in the limit T K — » 0, the form of 
P{T K ) is universal in the insulating phase. We find that 
P{T K ) ~ A T,7 Q where A is a non-universal amplitude and 
the exponent a depends only on the dimensionality. Fur- 
thermore, the same behavior is obtained for the metallic 
phase close to the metal insulator transition, with an iden- 
tical exponent. 

We discuss the problem of a single magnetic impurity 
Kondo-coupled to a disordered bath of non- interacting elec- 



trons. We represent the host by an Anderson model with 
random site energies uniformly distributed in an interval 
{—W/2,W/2}. Without loss of generality we shall set the 
chemical potential to zero and place the impurity at R=0. 
We shall denote the dimensionality of the system by D ( D 
= 1,2 and 3). The host is an insulator for W > W c where 
W c is finite in three dimensions (3D) and W c — in one 
(ID) and two dimensions (2D). 

Following previous authors we use slave-boson mean 
field theory to write down an implicit equation for T K 
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tanh(e/2T K ) p(e,0) 
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Here, p(e,0) = £^(s - e t ) \M°)\ 2 is the LDOS of the 
host at R = 0, where il>e(R) is the £-th eigenfunction of the 
host and eg is its associated eigenvalue for a fixed realiza- 
tion of the disorder. This formula is identical to what is 
obtained defining T K as the temperature at which second 
order perturbation in J breaks down [l(Tj |. 

In the first part of the Letter we discuss extensively the 
properties of P(T K ) in the Anderson insulator for which 
a number of results can be obtained using physical argu- 
ments and the solving analytically a simplified model that 
captures the essential physics. Then, we present the results 
of numerical simulations for all the insulating phases and 
for the metallic phase in 3D. 

In an Anderson insulator the wave functions are localized 
around random positions R^. At long distances from their 
peak positions they decay as |^(R — R^ )| 2 ~ exp(— 2|R — 
R-fl/Cf) where ^ = £(ee) is the energy-dependent localiza- 
tion length. Since £ (e) is smooth on the scale of T K we shall 
neglect its energy dependence and denote its value at e = 
by £. The eigenvalues Si are random variables and their 
distribution is Poissonian (llj . 

For a fixed configuration of the site energies the LDOS 
consists of a dense set of ^-functions with varying ampli- 
tudes. There are 0(£ D ) peaks with weights 0(^ D ) that 
come from states localized within a volume of size £ D 
around the impurity and infinitely many small peaks with 
exponentially smaller weights that come from states local- 
ized outside this volume. The mean separation between 
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the main peaks is AE^ ~ B(W)/£ D where B(W) is the 
bandwidth (i.e. the support of the LDOS). 

Consider first the contribution of the largest peaks to 
the LDOS and neglect the rest. The conditions of valid- 
ity of this approximation will be discussed below. It is 
convenient to split the interval of integration in Eq. (|TJ in 
subintervals of width T K << B(W). Within each of these 
the function tanh(e/2T K )/e is essentially constant and can 
be safely taken out of the sub-integral. Then, the LDOS 
in Eq. Q can be well approximated by the coarsegrained 
function p s {e) ~ N p (e,T K )/(£ D T K ) where N p (e, T K ) is the 
number of large peaks with energies in an interval of width 
T K around e. It is clear that N p {e,T K ) = £ D T K /B(W) 
and ~p s = B~ 1 (W). The LDOS may be split in two parts, 
p s (e) = p s + 5p s (e) where Sp s (e) is a random function with 
zero mean. 

Suppose that N p ^> 1. Then, 5p s (e) is a Gaussian ran- 
dom function with Sp s (e) Sp s (e') — p s ^ D S(e — e'). 
Equation Q and elementary properties of Gaussian ran- 
dom functions lead readily to 



P{T K ) oc exp 
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This is a log-normal distribution centered at T°, the Kondo 
temperature of a clean system with a uniform LDOS ~p s . 
The variance of the distribution is a 2 — (£k/£) where 
£k = l/(T°po) is the volume of the Kondo screening 
cloud H3. Since cr" 2 ~ TVp, Eq. © is valid for £ > £ K . 
In this regime most of the impurities have Kondo temper- 
atures in a narrow interval around T® and yet we are in an 
insulating state. 

This is not the end of the story. Even for £ > £ K the 
condition Np~ » 1 is violated for T K < B(W)/ £ D . For 
these temperatures, the number of large peaks with energies 
in an interval T K around the Fermi energy is small. They 
no longer dominate the LDOS as assumed and screening 
of the impurity must be done by electron states localized 
far away from the former. But these are precisely those 
that we neglected from the start. We then expect that for 
£ > £ K P{T k ) will show, in addition of the peak around T°, 
a low temperature tail appearing below T K ~ B(W)/ £ D 
carrying a small fraction of the weight of the former. For 
£ < £ K there is no temperature range in which Eq. @ holds 
and we do not expect to find a sharp peak near T°. The 
tail should still appear below ~ T°. 

In order to gain insights on the properties of the low- 
temperature tail in P(T K ) we now discuss an exactly solv- 
able "toy" model for the insulating phase that reveals in- 
teresting features that we shall retrieve in the numerical 
solution of the Anderson model presented below. 

As stated above the localized functions i/je(R) decay ex- 
ponentially away from their peak positions. However, they 
fluctuate randomly around the envelope ~ exp(— R/£)- In 
our toy model these fluctuations are ignored and the true 
eigenfunctions are replaced by normalized exponentials cen- 



tered around their random locations. Furthermore, these 
positions are assumed to be uncorrelated with the ener- 
gies. This holds deep in the insulating phase. Equation @J 
becomes 
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where the normalization Md = 2^ D fln£, D (D — 1)!, VLd is 
the solid angle in D dimensions and the last term in the 
above equation defines the random function Ft- 

A key feature of Eq. JJjJ is that Ft is a monotonically de- 
creasing function of T. A moment's reflection will convince 
the reader that this implies 



P{F T < 2/ J) = P(T K <T) + P{T K = 0) 



(4) 



where the first term on the right hand side is the integrated 
probability of T K from zero to T, and the second is the 
probability of the impurity remaining unscreened at T — 0. 
Therefore, the distribution of non-zero Kondo temperatures 
is given by 



P(T K ) 



2/J dW T {x) 



dT 



dx, 



(5) 



T=T K 



where Wt(x) — 5 (Ft ({si}) — x) is the probability density 
of Ft- In our toy model Wt(x) can be calculated analyti- 
cally leading to a closed form expression for P(T K ). 

Here, we are interested in the asymptotic properties 
in three relevant regimes. For T K /W l/Af-p ~ 
we find that P{T K ) sharply decays for large T K as ~ 
exp(— C£ D T K In 2 T K ) where C is a constant of order one. 
The same expansion is valid for £ > £ K , T K ~ T° and in 
this case we retrieve precisely Eq. (J5J. 

More importantly, we can now have access to the low- 
temperature region that was beyond the scope of the quali- 
tative analysis presented above. In the limit T — > we find 



P(Tk) 



|log(T K /H0r for T K ^0, (6) 



with a = D — 1. The asymptotic form of P(T K ) for this 
model is thus a power of the logarithm with an exponent 
depending only on dimensionality. The prefactor, propor- 
tional to £ D , decreases with increasing disorder. A more 
detailed analysis shows that this divergent contribution ap- 
pears at very low temperatures for £ 3> £k- With increasing 
disorder the tail shifts to higher temperatures. 

The model discussed above is based on a crude represen- 
tation of the localized states. Furthermore, it can obviously 
not describe the metallic phase in 3D or the metal-insulator 
transition (MIT). Therefore, we turn to our numerical re- 
sults for the Anderson model (AM) . It will be seen that the 
two models share the same qualitative properties. Namely, 
a divergent P(T K ) at low T K with the nature of the di- 
vergence depending only the spatial dimension and not on 
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disorder. We shall also see that P(T K ) acquires the same 
form in the localized and metallic phases in 3D. 

Our numerical calculations were performed using a rep- 
resentation of the host by the Anderson model on a regular 
lattice of size L D for D = 1, 2, and 3: 

H A = £ Ao C io ~ t 4a C ja > ( 7 ) 

where £j are random variables and the hopping integral 
between nearest neighbors t = 1 is the unit of energy. We 
show results for J = 1/4B, where B is the bandwidth of 
the clean system, and for Si uniformly distributed in the 
interval [—W/2,W/2]. We have checked the robustness of 
our results, in particular the value of the exponent of the 
power law, analyzing other values of J and considering a 
Gaussian distribution for Ej. 

For each lattice size and each realization of the random 
potential the Anderson Hamiltonian was diagonalized and 
Eq. JQ| was solved for T K using the numerical eigenfunctions 
and eigenvalues. Histograms were constructed repeating 
the procedure for many realizations of the potential and 
positions of the impurity. The size of the lattice and the 
number of realizations was increased until convergence of 
the histograms. The values of the parameters appear in the 
captions of the figures. 
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FIG. 1: P (T K ) in ID for several values of the amplitude of the 
disorder. Results for a chain of N s = 1000 sites. The num- 
ber of realizations of the potential is 1000. Inset: Fraction of 
unscreened impurities at T — as a function of disorder. 

Figure^displays the distribution of Kondo temperatures 
in D = 1. It is seen that at small disorder P(T K ) is strongly 
peaked at T°. For W -C 1 the shape of the peak fits accu- 
rately Eq. (j2J. With increasing disorder P(T K ) first broad- 
ens and departs from the log-normal distribution, and a 
tail appears at low temperatures. Upon increasing further 
the disorder, the peak at T K = disappears and P(T K ) 
becomes flat as in the toy model. 

The fraction of unscreened impurities at T = as a func- 
tion of disorder is shown in the inset to Fig. ^ P(T K = 0) 
increases continuously with disorder from vanishingly small 




FIG. 2: P(T K ) in 2D for several values of the amplitude of the 
disorder. Results for square lattices of -/V s = 2025 sites and 
1000 realizations of the amplitude of the potential. The slope 
of the thick straight line is -0.5. Inset: Fraction of unscreened 
impurities at T = as a function of disorder. 

values for W ~ to 0(1) for large W. We independently 
estimated £ as a function of disorder and checked that the 
crossover between the two regimes takes place when £ <~ £ K , 
as expected from our qualitative arguments. 

The results for the two dimensional case are displayed in 
Fig. El The same general trends are also observed in this 
case but in two dimensions P(T K ) is no longer flat but di- 
vergent at low T K as expected from our analytical results. 
Note that the straight lines corresponding to different val- 
ues of W are parallel meaning that P{T K ) is a power law 
with a disorder independent exponent. The thick straight 
line has slope —1/2 suggesting that a — 1/2 in 2D. 

For temperatures below the range plotted in the figure 
finite-size effects come into play. These are expected for 
T K < A, the typical level spacing at the Fermi level oc 
L~ D 13]. We have thus performed a careful study of these 
effects and checked that the power law behavior extends 
to lower and lower values of T K as the size of the system 
increases. 

The results for the 3D case, displayed in Fig. are the 
most interesting ones because in this case there is a MIT 
at W c ~ 16.5 14]. As in the 2D case, P(T K ) is a power 
law at low temperature but the slope of the thick straight 
line is now —3/4 suggesting that this is the value of the 
exponent in 3D. A striking new feature in 3D is that the 
form of P (T K ) and the exponent are the same on both sides 
of the MIT. This can be seen by comparing the two low- 
ermost curves in Fig. that correspond to W > W c , to 
the two above them that correspond to W < W c . The to- 
tal contribution of the impurities to the susceptibility or 
specific heat can be calculated using interpolation formulas 
(see e.g. Ref. 0) and weighting the individual contribu- 
tions with P(T#). An important consequence is that the 
magnetic susceptibility and the specific heat coefficient in 
this model should exhibit a power law behavior with a dis- 
order independent exponent, x ~ 7 ~ T~ a , and be 
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FIG. 3: P(T K ) in 3D for several values of the amplitude of the 
disorder. Results for cubic lattices of N s — 2197 sites and 1000 
realizations of the amplitude of the potential. The MIT occurs 
at W c ~ 16.5. The slope of the thick straight line is -0.75. Inset: 
size and disorder dependence of P(T K — 0) in the metallic phase. 



cess in that context A power law P(T K ) was found 

for the disordered Anderson model on a Bethe lattice us- 
ing DMFT in Ref. An important difference between 
this work and ours is that the exponent of the power law 
found in Ref. depends on the strength of disorder and 
is therefore non-universal. It would be very interesting to 
investigate how our results for the diluted case evolve with 
increasing density and what is the role of the lattice geom- 
etry. 
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After the completion of this work we became aware of a 
related study by S. Kettemann and E. R. Mucciolo 



continuous across the MIT. Interestingly, this is precisely 
what was observed in Si:P which has a MIT as a function of 
the concentration of P and shows signatures of local mag- 
netic moments in both the metallic and insulating phases 
0,0]. Although the agreement obtained is very appealing, 
this system contains physical ingredients which are out of 
the scope of our model. 

We found that in our finite systems a fraction of the spins 
remains unscreened down to zero temperature. This frac- 
tion, very small for small disorder, increases as W — > W c . 
Since on general grounds, in the thermodynamic limit, all 
spins must be screened at T = in the metallic phase, the 
observed effect is due to finite-size. The inset to Fig. dis- 
plays the results of a finite-size scaling analysis of P(T K = 0) 
for several values of W. It can be seen that P(T K = 0) 
scales to zero in the metallic phase with increasing L as 
it should. Moreover, P(T K = 0) oc L~@( w ' with an expo- 
nent (3(W) that depends on the disorder and that decreases 
as W increases. In contrast, deep in the insulating phase, 
P(T K = 0) remains finite in the thermodynamic limit Q. 

In this Letter we have studied the single-impurity prob- 
lem and our results are strictly applicable in the diluted 
limit. We showed that the structure of the LDOS induces 
a power law P(T K ) in the insulating phase. This low T K 
divergence is associated to impurities screened by states lo- 
calized at distances larger the localization length f. The 
statistical properties of these far states are such that a 
change in the disorder strength shifts the onset temper- 
ature of the diverging behavior but not the nature of the 
divergence. Surprisingly, the exponent of the power law is 
disorder independent and in 3D the same behavior persists 
across the MIT with an identical exponent. 

In concentrated systems interaction effects can not be ig- 
nored. An extension of the Dynamical Mean Field Theory 
(DMFT) method for random systems was used with suc- 
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